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Abstract 

We test the internal consistencies and the reliability of the existing estimates of the decay 
constant fs in the static limit, the meson-quark mass gap A and the kinetic energy K 
of a heavy quark obtained from the heavy quark effective theory (HQET)-sum rules. 
Finite energy local duality sum rules (FESR) have also been used to fix approximatively 
the value of the continuum energy and to study the correlations among these different 
parameters. Then, we deduce to two-loop accuracy: A = (0.65 ± 0.05) GeV, K = 
-(0.5 ± 0.2) GeV 2 , implying the value of the pole mass in HQET: M b = (4.61 ± 0.05) 
GeV. By combining the results from the sum rules in HQET and in the full theory, we 
obtain /J? = (1.98 ± 0.31)/^ and the quark mass dependence of the pseudoscalar decay 
constant: f P VW = (0.33 ± 0.06) GeV 3 / 2 ^ 1 {l - 2a s /3vr - 1.1/Mq + 0.7/M 2 } . 
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There has been a lot of effort during these last two years for estimating the parameters 
of the heavy quark effective theory (HQET) [f] using QCD spectral sum rules (QSSR) 
[2, 3] in the infinite mass limit [4]-[14]. Among other existing estimates, we shall focus 
our discussion on the determinations of the decay constant fs in the static limit, the 
meson-quark mass gap A and the kinetic energy K of a heavy quark. Our aim is to test 
the reliability of the existing estimates of these quantities and to try to extract the most 
reliable information from the sum rules analysis. In so doing, we shall first reconsider 
the HQET sum rules analysis by using an update of the QCD parameters and by testing 
the stability of the results with the variation of the Laplace sum rule variable r and/or 
of the continuum threshold E c 1 . If this stability exists, then, we shall compare it with 
the corresponding conservative values of the set (r, E c ) obtained from QCD spectral 
sum rules (QSSR) in the full theory 2 . We shall also use local duality finite energy sum 
rule (FESR) to fix approximatively E c and study the different correlations among these 
different parameters. 



2 Definitions and normalizations 

The decay constant of the pseudoscalar meson is defined as: 

{0\J 5 "\B) = V2f B q" A = 93.3MeV, (I) 
where J$ is the axial current of the full QCD theory. In HQET, the decay constant reads: 

m\B) = -j=f B v», (2) 

to which corresponds the RGI coupling: 

f B =Mv)(a s f^(l-S^y (3) 

where f3i = —1/2(11 — 2nj/3) is the first coefficient of the /3-function and 6 ~ —0.23 for 
nj = 4 flavours; v^ is the heavy quark velocity; J$ is the axial current in HQET, which 
is related to the former J$ of the full theory by: 



1 The meaning of this stability has been discussed in [9] in quantum mechanics , by taking the Laplace 
transform of the eigenvalue solution of an harmonic oscillator and by comparing the approximate series 
with the known complete solution. The optimal result is obtained at the minimum or at the inflexion 
point of the curves. 

2 These conservative values correspond to the value of E c comprised between that at which the r- 
stability is reached and that where we have ^-stability. These criteria include all different alternative 
criteria used in the literature and gives a more precise meaning of the so-called sum rule window of SVZ. 
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The meson-quark mass gap A is defined as [1] 3 : 

Mb = M 6 + A-A-(^ + 3E), (6) 

where: 

K = ^-(B(v)\lC\B(v)) and E = -L-( B (v)\S\B(v)) (7) 

corresponds respectively to the matrix elements of the kinetic and chromomagnetic oper- 
ators: 

K = h{iDfh and S = hia^F^h, (8) 

where h is the heavy quark field and F^ v the electric field tensor. 



3 HQET sum rules 

Let us start our discussion from the relativistic Laplace sum rule (LSR) [2, 3]: 

Cr{t t ) = dt e~ tTr -ImII(t), (9) 

JO 7T 

where t C} t t are respectively the QCD continuum threshold and the Laplace sum rule 
variable, which a priori are free external parameters in the sum rule analysis; LI is a 
generic notation of the hadronic Green's function. HQET sum rules correspond to the 
limit of LSR where the heavy quark mass tends to infinity. In this case, it is more 
convenient to introduce the non-relativistic variables r and E, which are related to the 
former as: 

t = (E + M Q )\ r e = r r ^ (l + O (^-)) , (10) 

or: 

uj = 2E } T^ = T r M Q} (11) 

which remain finite when Mq tends to infinity. These non-relativistic sum rules have been 
discussed earlier [9, 10], and have been exploited in the context of HQET in [4]-[8]. They 
have also been used in full QCD but for increasing the value of Mq in [1 1]— [14] in order 
to reach the HQET regime from below. It has also been shown [12] that the Laplace and 
moments sum rules are equivalent, when the number of derivatives n of the n-th moments 
and Mq are large, but their ratio r = u/Mq remains finite 4 . Therefore, without loss of 



3 We are aware of the fact that in the lattice calculations, A defined in this way can be affected by 
renormalons [15]. 

4 However, it has to be noticed that in the analysis of the heavy-to-light transitions, the optimal results 
are obtained at small n- values [17], such that, the Moments sum rule does not tend to the Laplace one, 



when Mq is large. 
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rule: 



or equivalently with: 



£(t 6 ) = dE e~ ETe —lmIl(E) } (12) 

7T 



f Wc 1 

C(t u )= dto e-^-Imn(cu); (13) 

JO 7T 



both sum rules have been used in the literature. The apparent advantage of Eq. (13) 
is that it optimizes at a value of t w ~ T e /2, which is a trivial fact due to their intrisic 
definitions, but it is psychologically more comfortable to reach the stability at a value of 
t -1 around 1 but not 2 GeV -1 , therefore giving a better justification for the convergence 
of the Operator Product Expansion. 



4 The mass gap A and the decay constant /# 

The sum rule for the decay constant and for the mass gap A is well-known in the literature. 
It reads [4]-[6], [3, 12]: 

-™4 + ^( j+ 4H^'}- 

(14) 

where r = t w ; A = 704/9 — 1127r 2 /27 [5], but its effect is not significant in the analysis. 
The subtraction point v can be naturally chosen at v = 1/r, as dictated by the RGE 
which governs the Laplace sum rule [18]. However, in order to study the //-dependence 
of the results, we shall take it from 1 GeV to 2/r e , which covers the choice done in the 
literature. The mass gap can be deduced from the log-derivative of Eq. (14) with respect 
to r. We shall use the values of the QCD scale for 4 flavours to two-loops: 

A 4 = (317 ± 84) MeV, (15) 

corresponding to a s (Mz) — 0.118 ± 0.006 [19], and of the QCD non-perturbative param- 
eters [3]: 

Ml = (0.8±0.01) GeV 2 {uu)(l GeV) ~ -(0.223 GeV) 3 {a s G 2 ) = (0.06±0.03) GeV 4 . 

(16) 

We study the r-behaviours of fs and of A in Figs, la, c for a fixed value of E c . At the 
approximation where the correlator is computed, there is a good stability for the values 
of r in the range of 1 to 2 GeV -1 . The dependence of these optima on E c are shown in 
Fig. lb, d, where fs has a much better stability than A. Taking the conservative range 
of E c suggested by the analysis of fs from QSSR in the full theory [3], [1 1]— [14] : 

E c ~ 1 ~ 1.6 GeV, (17) 
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f B ~ (0.21 ~ 0.28) GeV 3/2 A ~ (0.52 ~ 0.70) GeV. 



(18) 



The weak dependence of the results on E c shows the irrelevance of the criticisms of [7] on 
the choice of E c done in [11], while the ones of [16] do not concern the working regions of 
the sum rules. It should be noticed that the effect of the subtraction scale in the range 
proposed previously is not important. Varying a s within the error bar given previously, 
one affects A and fs respectively by 2 and 7%. The error due to the quark condensate is 
much smaller. The most significant error not taken into account in the HQET-literature 
with the exception of [12] is the one due to the unknown a 2 -term. This is important for 
fs as the perturbative radiative correction is huge, which is of the order of the lowest- 
order term in the two-point correlator. In the absence of an evaluation of this term, a 
reasonable estimate is obtained by assuming that its coefficient grows like a geometric 
sum or, more conservatively that its effect cannot be as large as the known a s -term. In 
this way, an error of about 25 % on fs is introduced 5 . Then: 

f B ~ (0.25 ± 0.04 ± 0.06) GeV 3/2 , (19) 

where the first error is due to the sum rule procedure and the second one is mainly 
dominated by the unknown a 2 -term. The values of fs and of A are in agreement with 
the previous results in [5, 6], but the error for fs must have been underestimated there. 
We show in Fig. Id the result of [6] (dashed line). The small difference comes from the 
different value of A 4 used in the analysis. Moreover, these values are also in agreement 
with the ones from the sum rule in the full theory, obtained by increasing the quark 
mass-value [12]: 

f B ~ (0.33 ± 0.06) GeV 3/2 , A ~ (0.6 ~ 0.8) GeV, (20) 

where fs corresponds to: 

(Jb\[Mb^° ^ (0.42 ± 0.07) GeV 3/2 . (21) 



5 The kinetic energy K 

The sum rules for the kinetic energy have been proposed in [6, 8], where the two-point 
and three-point functions have been used respectively. However, although the two-point 
function is a convenient quantity to work with, one cannot disentangle there the different 
contributions, which are of the same order in 1/Mq. Therefore, we shall only consider 
the three-point function analysis from [8]. The sum rule reads: 



f 2 B Ke-^ ~ £ C dE J^dE' e -( E+E >' 2 U-^E 4 6(E-E') 



aJv) fil 4tt 2 , ^ N 
1 + — + — 2 In Et 



TT 
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3 In the following this error will be systematically included. 
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— — — ^ + u ^- + u- 

7T 7T J 
+ 4o £ H r + + + ^ _ | 

3 7T L J 

-L^Q^E _ ^ + | aa -2/(3ft) M 2( Su )( I/ ) 5 ( E ) 5 ( E ') 

47T 8 

(22) 

where r = r e . The kinetic energy K can be obtained by taking the ratio of the sum 
rules in Eqs. (22) and (14). A similar analysis can be done here, which is summarized in 
Fig. 2a, b. Again a good stability in r is reached for E c > 0.8 GeV. This suggests the 
phenomenological but not rigorous lower bound: 

\K\ > 0.4 GeV 2 . (23) 

However, the dependence on E c is violent. Hopefully, a slight inflexion point appears for 
E c ~ 1 ~ 1.2 GeV, at which one can extract the optimal result: 

K ~ -(0.5 ~ 0.6) GeV 2 , (24) 

but our confidence level in this estimate is less than for the previous parameters. A more 
conservative attitude is to consider the whole range spanned in Eqs. (23) and (24) and 
to enlarge the error by a factor 2. Thus, one obtains: 

K ~ -(0.5 ± 0.2) GeV 2 . (25) 

Though in agreement with that of [8], this result suggests that a lower value of K about 
-0.3 GeV 2 is not absolutely excluded. 



6 Local duality FESR constraints 

Let us now test the consistency of the previous estimates by using local duality FESR. 
The constraints from the lowest FESR moments can be derived from the LSR by taking 
the limit r — > 0. In this way, one can deduce the set of approximate equations: 

A « ^1 + £!<««) j 

- - 14-T^Hf)}- 

We solve these constraints by using as input the value of A obtained previously, which 
is quite stable versus the variation of E c and which is also less affected by the radiative 
correction than This gives: 

E c w 1 GeV (27) 
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] B « (0.40 ± 0.07) GeV 3/2 
K « -0.41 GeV 2 , 



(28) 



which agrees within the errors with the previous estimates and indicates the consis- 
tency of the whole results from the LSR and FESR. 

7 Conclusions 

Using the LSR and local duality FESR, we have derived the values of the HQET param- 
eters, at which the results are consistent with each other. We conclude that, besides the 
eventual contribution of renormalons, the most reliable estimate is that of A, although it 
is slightly affected by the variation of E c . Indeed, contrary to its radiative correction 
is quite small (ratio of sum rules). We can consider as the best estimate, the interaction 
region between the HQET and full QCD sum rules results in Eqs. (18) and (20), which 
suggests: 

A ~ (0.65 ± 0.05) GeV. (29) 

This value is slightly higher than the existing estimates [6, 5] and the recent average [22], 
which is based on a more precise value A ~ (0.55 ± 0.05) GeV from HQET sum rule 
than the one in Eq. (18). For the kinetic energy, the result is inaccurate due to the huge 
dependence of the estimate on the change of the value of E c . Our conservative estimate 
in Eq. (25), which spans the range of values provided by the LSR and FESR methods is 

6. 

K ~ -(0.5 ± 0.2) GeV 2 . (30) 

By combining the previous results with the well-known estimate of the chromomagnetic 
matrix element to leading order in 1 / M b : 

E ~ i(M 2 « - M 2 B ) = 0.119 GeV 2 , (31) 

one can deduce from Eq. (6), the value of the quark mass to two-loop accuracy entering 
in HQET: 

M b ~ (4.61 ± 0.05) GeV. (32) 

This value is in agreement (within the errors), with previous analysis from HQET-sum 
rules [6] and with the most recent determinations of the pole mass to two loops from the 
relativistic sum rules fit of the B and B* masses in the full theory [22]: 

M b ~ (4.63 ± 0.08) GeV, (33) 

and from a non-relativistic sum rule analysis of the bottomium systems: 

M b ~ (4.69 ± 0.03) GeV, (34) 

6 We are aware of the importance of this quantity in the analysis of the inclusive B-meson decays [21]. 
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contributions in 1/M& to the relation given in Eq. (6). It may also indicate that the still 
quantitatively less controlled renormalon effects to the quark mass defined in HQET[23], 
is small regardless the definition of A introduced in the analysis. 7 . For the decay constant 
/b, the radiative correction is huge, and it amounts to about 100% of the lowest-order 
contribution in the corresponding two-point correlator. Therefore, in the absence of an 
evaluation of the a 2 -term, one can only suggest the conservative range spanned by the 
results of the two different HQET sum rules. Then, we obtain 

fp ~ (0.34 ±0.16) GeV 3/2 , (35) 

where, as already mentioned, the error is mainly due to the estimate of the unknown 
a 2 -corrections. This result can be transformed into the value of the decay constant in the 
so-called static limit, by using the relation in Eq. (5). Therefore, one obtains: 

(2.02 ±0.9)./;. (36) 

This result is in agreement with the value: 

(1.97 ±0.33)./; (37) 

obtained from the sum rule in the full theory [12], which does not however suffer from a 
large a s -correction and therefore has a much smaller error. One can combine these two 
results and quote the final value of the decay constant in the static limit from the sum 
rules: 

fp° ~ (1.98 ± 0.31)./; and f B ~ (0.33 ± 0.06) GeV 3/2 . (38) 

This result is in good agreement with the present lattice results /J? ~ (1.74 ± 0.27)./; 
compiled by [20]. This result can be used for fixing the coefficients A and B of the 1/Mq 
and 1/Mq corrections on fp from Eq. (5): 

W^^.(M,))^(l-f^-A + .*), m 

In so doing, we shall use the values [22, 12, 24]: 

fp> — (1.37 ± 0.04 ± 0.06)./; and fp — (1.49 ± 0.06 ± 0.05)./;. (40) 

First, we use a linear fit and fix simply A from the ratio fp/ fp. We obtain: 

A= (0.6 ±0.1) GeV, (41) 

which agrees within the errors with the previous findings [11], [5]-[8], [12]. Secondly, we 
determine A and B by solving numerically these three equations. We obtain: 

A « 1.1 GeV and B « 0.7 GeV 2 , (42) 

where A and B have large errors. These two terms agree within the errors with the earlier 
results [11], which suggests that the 1/Mq correction is smaller than the 1/Mq one even 
at the charm quark mass, been 



7 An indication of the smallness of this effect can also be inferred from the comparison of the values of 
the pole mass from relativistic and non-relativistic bottomium sum rules [22] 



7 



Fig. 1 : Analysis of the r- (a, c)) and E c - (b, d) behaviours of fs and A for A 4 = 0.317 
GeV. The dashed line is the prediction of [6] using A 4 = 0.25 GeV. Fig. 2 : The same 

as Fig. 1 but for K. In a), we show the effect of the subtraction point. 
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